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Abstract 

The usual procedure of including a finite number of vertices in Non Perturbative Renormalization 
Group equations in order to obtain n-point correlation functions at finite momenta is analyzed. 
This is done by exploiting a general method recently introduced which includes simultaneously all 
vertices although approximating their momentum dependence. The study is performed using the 
self-energy of the tridimensional scalar model at criticality. At least in this example, low order 
truncations miss quantities as the critical exponent n by as much as 60%. However, if one goes to 
high order truncations the procedure seems to converge rapidly. 
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I. INTRODUCTION 



In nearly all fields in physics, there are systems having a large number of strongly cor- 
related constituents. These cannot be treated with usual perturbative methods. Phase 
transitions and critical phenomena, disordered systems, strongly correlated electrons, quan- 
tum chromo dynamics at large distances, are just a few examples which demand a general 
and efficient method to treat non-perturbative situations. In problems as those just quoted, 
the calculation of correlation functions of the configuration variables is, in general, a very 
complicated task. 

The non perturbative renormalization group (NPRG) has proven to be a 

powerful tool to achieve this goal. It presents itself as an infinite hierarchy of flow equations 
relating sequentially the various n-point functions. It has been successfully applied in many 
differentproblems, either in condensed matter, particle or nuclear physics (for reviews, see 
e.g. 0, Q, Q]; a pedagogical introduction can be found in Q]). In most of these problems 
however, one is interested in observables dominated by long wavelength modes. In these 
cases, it is then possible to approximately close the infinite hierarchy of NPRG equations 
performing an expansion in the number of derivatives of the field. This approximation 
scheme is known as the derivative expansion (DE) 10|. The price to pay is that the n-point 
functions can be calculated only at small external momenta, i.e. smaller than the smallest 
mass on the problem (vanishing momenta in the case of critical phenomena). 

In many other physical problems however, this is not enough: the full knowledge of the 
momentum dependence of correlation functions is needed in order to calculate quantities of 
physical interest (e.g. to get the spectrum of excitations, the shape of a Fermi surface, the 
scattering matrix, etc.). There have been many attempts to solve the infinite system of flow 



equations at finite momenta; most of them are basec 



by Weinberg 



111 ]. Although some of these attempts 



on various forms of an early proposal 
Q introduce sophisticated 
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ansatz for the unknown correlation functions appearing in a given flow equation, most efforts 
simply ignore high order vertices. In all these works, only low order vertices are taken into 
account: usual calculations do not even include the complete flow of the 3- and 4-point 
functions. Moreover, it is not possible a priori to gauge the quality of such approximations 
schemes. 

Recently, an alternative general method to get n-point functions at any finite momenta 
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within the NPRG has been proposed It has many similarities with DE. First, it is 

an approximation scheme that can be systematically improved. Second, the scheme yields 
a closed set of flow equations including simultaneously an infinite number of vertices; one 
thus goes far beyond schemes including a small number of vertices, as those quoted in 
the previous paragraph. Moreover, it has been proven 16] that in their corresponding 
limits, both perturvative and DE results are recovered; this remains valid at each order of 
the respective expansion. Finally, in the large- N limit of O(N) models, the leading-order 
(LO) of the approximation scheme becomes exact for all n-point functions. (The expression 
"leading order" means the first step in the approximation scheme; it does not refer to an 
expansion in a small parameter which usually does not exist in these kind of problems.). 



In 



171 ]. the method has been applied, in its leading order, to the calculation of the self- 



energy of the scalar model, at criticality. That is, we have fine-tuned the bare mass of the 
model in order for the correlation length to be infinite, and then we have studied the full 
range of momenta, from the high momenta Gaussian regime to the low momenta scaling 
one. At this order of the approximation scheme the self-energy is expected to include all 
one loop contributions and to achieve DE at next-to-leading order (NLO) precision, in the 
corresponding limit [if]]. The numerical solution found in [l^] verifies these properties. 
Moreover, the function has the expected physical properties in all momenta regime. First, 
it presents the correct scaling behavior in the infrared limit. The model reproduces critical 
exponent rj with a level of precision comparable to the DE at NLO. Moreover, contrarily to 
DE, the anomalous power-law behavior can be read directly from the momentum dependence 
of the 2-point function. Second, it shows the expected logarithmic shape of the perturbative 
regime even though the coefficient in front of the logarithm, which is a 2-loops quantity, 
is only reproduced with an error of 8%. In order to check the quality of the solution in 
the intermediate momentum region, a quantity sensitive to this crossover sector has been 
calculated: one gets a result almost within the error bars of both Monte-Carlo and resumed 
7-loops perturbative calculations. Please observe that this quantity is extremely difficult to 
calculate: even these sophisticated methods give an error of the order of 10%. 



Another interesting similarity between DE and the method presented in 16] is that, as a 
price to pay in order to close the equations including an infinite number of vertices, one has 
to study the problem in an external constant field. Accordingly, one ends up with partial 
differential equations which may be difficult to solve. A useful approximation scheme, widely 
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used in DE calculations, is to perform, on top of the expansion in derivatives of the field, 
an extra expansion in powers of the field (see e.g. [3]), in the spirit of Weinberg proposal. 



During the last 10 years, this strategy has been widely used 
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in many studied situations this expansion seems to converge (generally oscillating 
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23l . [2J, |25|, while in many others it does not 22j, |26(. In d = 2, the field expansion has 



2d, 



been explored with no indication of convergence for critical exponents, even going to high 



orders 



23]. 



In this work we shall explore this procedure of expansion in powers of the field, in the 
framework of the calculation scheme presented in 16| . More precisely, we shall make a field 
expansion on top of the already approximated 2-point function flow equation solved in 1^] . 
Then we shall compare results with and without field expansion. In doing so, we have two 
goals. First, we shall study the apparent convergence of this procedure. This comparison is 
essential if one hopes to apply the scheme described in [3] to situations more complicated 
than that considered in 1_7|] . For example, within DE scheme, when trying to go to higher 
orders or when considering more involved models, the expansion in powers of the field on top 
of the corresponding approximate flow equations is sometimes the only practical strategy to 



solve them 



23 



28] . The second and more important goal is the following: as we shall see 



in section 1111} truncation in powers of the field is equivalent to ignoring high order vertices 
in the flow equations. Thus, the comparison presented here can help to estimate the quality 
of the calculations made so far to get n-point functions at finite momenta neglecting high 
order vertices. 

The article is organized as follows. In the next section we describe the basics ingredients 
of both the NPRG and the approximation scheme introduced in [3] . We also present the 



results obtained in 17|, when this scheme is used to find the 2-point function of the scalar 



model. In section IIIH we apply the expansion in the field at various orders and compare 
these results with those found in [ljj]. Finally, we present the conclusions of the study. 



II. GENERAL CONSIDERATIONS 



Let us consider a scalar field theory with the classical action 

1 



S 



• J ^(^)) 2 + ^ 2 (x) + |^(x) 



d a x 



(1) 
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Here, r and u are the microscopic mass and coupling, respectively. 

The NPRG builds a family of effective actions, T K [(j)] (where 4>(x) = (cp(x))j is the ex- 
pectation value of the field in presence of an external source J{x)), in which the magnitude 
of long wavelength fluctuations are controlled by an infrared regulator depending on a con- 
tinuous parameter k. One can write for T K [<j)] an exact flow equation j^, Q, 21, 29]: 

(2) 

where T K is the second functional derivative of T K with respect to 4>(x), and R K denotes 
a family of "cut-off functions" depending on k: R K {q) behaves like k 2 when q <C k and it 
vanishes rapidly when q ^> k Q, Q. The effective action T K [0] interpolates between the 
classical action obtained for k = A (where A -1 is the microscopic length scale), and the full 
effective action obtained when k — > 0, i.e., when all fluctuations are taken into account (see 
e.g. [31). 

By differentiating eq. (J2]) with respect to 4>(x), and then letting the field be constant, one 
gets the flow equation for the n-point function in a constant background field <p. For 
example, for the 2-point function one gets: 

d K Y^{p;<j>) = j -^ d d K R k {q)[G K {q^)^\v,q,-V-q-^) 
xG K (q + p; 4>)Yf{-p,p + q, -q; <p)G K {q; 0) 
-\G K (q; cf>)rW{p, -p, q, -q; 0)G K (g; 0)}, (3) 

where 

G-\q;<P)=rW(q,-q^) + R K (q 2 ), (4) 



and we used the definition 



(2n) d 6® f J>j r £°(Pi. ■■■iPnA) = J d d x, ... J d d x n e^U^ 



s n r k 



50(si) . ..5<j)(x r . 



(x)=4> 
(5) 



The flow equation for a given n-point function involves the n+1 and n + 2 point functions 
(see, e.g., eq. ([31)), so that the flow equations for all correlation functions constitute an infinite 
hierarchy of coupled equations. 



In 



161 ]. a general method to solve this infinite hierarchy was proposed. It exploits the 



smoothness of the regularized n-point functions, and the fact that the loop momentum q in 



5 



the right hand side of the flow equations (such as eq. (|2J) or eq. (j3j)) is limited to q ^ n due 
to the presence of d K R K (q). The leading order of the method presented in [3] thus consists 
in setting 

^ ( k\pi,P2, -,Pn-l + q,Pn - q) ~ T^ ) (p 1 ,p 2 , p„_i,p n ) (6) 

in the r.h.s. of the flow equations. After making this approximation, some momenta in some 
of the n-point functions vanish, and their expressions can then be obtained as derivatives of 
m-point functions (m < n) with respect to a constant background field. 

Specifically, in the flow equation for the 2-point function, eq. ([3]), after setting q = in 
the vertices of the r.h.s., the 3- and 4-point functions will contain one and two vanishing 
momenta, respectively. These can be related to the following derivatives of the 2-point 
function: 

r?(P,-P,o;t) = ar *\ p ^\ rw(p,- p ,o,o;0) = g^|f^ . (7) 

One then gets a closed equation for I« (p; <p): 

* r (2) f 2 A \ 7 (3)/ ,s (dT [2 \p, -p;(f>)\ 1 (2), ,x d 2 T (2 \p, -p;4>) , , 

Kd K r%>(p ;</>) = jy{p,K;4>) I — I --iy{K-,(f>) — 2 , (8) 

where 

J^{p-K-<j>) = J ^^ K ( g 2 )G (t (p + g ;0)G^- 1 )(g;0) ) (9) 



if \k; 0) = f A^ K d K R K (q 2 )G n K (q; 0). (10) 



and 

i d 
(27r) d 

In fact, in order to preserve the relation 

ri 2) ( P = O;0) = ^, (11) 

V R (<f>) = T K [<p(x) = 0]/Vol being the effective potential, it is better to make the approxima- 
tion (jSJ) (followed by (JTj)) in the flow equation for £ K (p; 0) defined as 

E K (p; 0) = rj?>(p; <P)-p 2 - rW(p = 0; 0). (12) 

The 2-point function is then obtained from T^ 2 '(p; <p) = d 2 V K (cf)) / d<p 2 + p 2 + S K (p; 0), which 
demands the simultaneous solution of the flow equations for V K (<p) and S K (p; </>). 
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As shown in 17|, even if the complete solution of these equations is a priori complicated, 



a simple, and still accurate, way of solving them consists in assuming in the various integrals 

G-\ q] 0) ~ Z K q 2 + d 2 VM/d<p 2 + R K (q 2 ), (13) 

where Z K = Z R (<fi — 0), with Z K (<j)) = l+<9£ K (p; 4>)/dp 2 \ p= o. This approximation is consistent 
with an improved version of the Local Potential Approximation (LPA, the first order of the 
DE), which includes explicitly a field renormalization factor Z K [7]. Doing so, the u p = 0" 
sector decouples from the p ^ one. Here, by "p = 0" we mean the sector describing vertices 
and derivative of vertices at zero momenta, i.e., flow equations for V K and Z K . Moreover, it 



is useful to use the regulator 



20] 



R K (q 2 ) = Z K (K 2 -q 2 )Q(K 2 -q 2 ), (14) 
which allows the functions J^\p', k; 4>) and 1^ (k; (j)) to be calculated analytically. The 



corresponding expressions can be found in 17J. In fact, all quantities are functions of 



p = 2 /2. The problem is then reduced to the solution of the three flow equations for V K (p) 
and Z K (p), for the p = sector, and for p), in the p ^ one. As only the "effective 



mass" 



<(p)s ^m = ?m +2p ?!m (15) 



(and its derivatives with respect to p) enters in the p ^ sector, it is more convenient to 
work with the flow equation for m 2 K (p) instead of that for V K (p) itself. The non-trivial fact is 
that by differentiating twice the flow equation for V K (p) w.r.t 0, one gets a closed equation 
for m 2 K {p). 

In order to make explicit the fixed point in the k — ► limit, it is necessary to work with 
dimensionless variables: 

fi K {p) = Z' 1 k- 2 m\{p) , xM = Z- 1 Z K (p) , p = K d l Z K K 2 ~ d p , (16) 

which, in the critical case, have a finite limit when k — > 0. Above, is a constant 
conveniently taken as K^ 1 = d 2 d ~ 1 K d l 2 F(d/2) (e.g., K 3 = l/(67r 2 )). In the p ^ sector, 
the dimensionful variable p in the self-energy flow equation makes S K (p; p) reach a finite value 
when k — > 0. As discussed in 1^], the inclusion of the flow equation for the renormalization 
factor Z K (p) is essential in order to preserve the correct scaling behavior of T( 2 \p; p) in the 



infrared limit. Doing so, in the critical case, the function T^ 2 \p; p)/[Z K K 2 ) has to reach a 
fixed point expression depending on p and p/tz, when K,p <?i and p ~ 1. 

Putting all together, in d — 3, the three flow equations that have to be solved are 

va fn' K (p) + 2pn'Z(p) WAP? 



Kd K p K (p) = -(2 - T] K )p K {P) + (1 + Vk)w'k(p) ~ ( 1 



5 J V (l + /i.(p)) 2 (l + /i K (p)) 3 



(17) 



and 



Kd K x K (p) = VkXk(p) + (1 + Vk)px' k (p) ~ 2 P (i+£ ( fj) ) 
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+ (i - f ) 8pxuP)n5fg)F - x fgry > 



together with 



X ^(0)/5 + (1 + ^(0))' 



^ = ,.,,nw c A - .. 77^2 > K<9 ^ = > ( 19 ) 



for the p = sector, and 



^( P ,p) = (i + Vk )pK(p,p) + 2 i0^ (wr P ) - $3$ 

(w«(p)K(p,p) + ^) - (S2# (K(P,P) + 2~pK(p,p)) (20) 

for the p one. In these equations, the prime means dp and we used the explicit expression 
for if = 2K 3 K 5 - 2n Z 1 K - n (l-ri re /5)/(l + /i K (p)) n . In eq. (120]) . we introduced the dimensionless 
expression f K defined as Jg (p; k; p) = K 3 /t _1 Z ft r 2 /(l + p K (p)) 2 x / re (p; p), with p = p/n. 

In 17] . this strategy is used to get the 2-point function of the scalar model at criticality 
and zero external field (i.e., S(p = 0, p = 0) = 0), in d = 3. As recalled above, the function 
thus obtained has the correct shape, either in the scaling, perturbative and intermediate 
momenta regimes. 



III. EXPANSION IN POWERS OF THE FIELD 

In this section, we shall compare the solution obtained in using the procedure de- 
scribed above, with the solution of the same three flow equations expanded in powers of 
p and truncating up to a given order. Before doing so, let us first consider only the flow 
equation for the potential or, equivalently, that for the effective mass, i.e., eq. (jT7|) . with 
Z K = 1 (r] K = 0). This corresponds to the pure LPA sector and it is thus independent of 

8 



the scheme presented in 16|. In d = 3, its expansion in powers o: 



' the field has been widely 



studied during the last ten years, using various regulators |22|, [24J . Recently, another inter- 
esting truncation scheme has also been considered in 30| showing much better convergence 
properties. However, here we shall consider the simpler expansion in powers of the fields; 
as shall be seen bellow this is the field expansion that can be compared to usual truncation 
in the number of vertices. It has been shown 25j that, using the regulator we consider here 



(see eq. {ED), this expansion seems to converge. This result follows when expanding both 
around finite and zero external field, although faster in the first case. In 25| the convergence 
in this situation has been discussed studying the critical exponent v. In order to strengthen 
this conclusion, as a first step in our study, we have analyzed the effect of the expansion on 
the function p K (p): 



oo ^ 



n) pn ^ 



(21) 



n=0 



More precisely, we shall gauge the impact of truncating this sum on the fixed point values of 
the coefficients pk , which are proportional to vertices at zero momenta and zero external 
field. This study is motivated by the fact that these pffl shall appear in the £ K (p;p) flow 
equation, eq. (1201) . when the later shall be expanded around p = 0. Results are shown 
in Figure [TJ The four plots present the fixed point value for the first 4 couplings, p^2 , 
n = 0, • • ■ ,3. For each coupling, we present the result which follows by solving the complete 
LPA equation, eq. ( JT7|) . together with the result obtained with the equation expanded in 
powers of p. For example, when going only up to the first order (i.e., neglecting all plif 1 with 
n > 2), the corresponding equations for p£ and Pk ^ , are: 

(1 - T)Jh)p^ 



(0) 



1 



/i 0) ) 2 



(22) 



and 



1 + /^ 3 



(23) 



which have to be solved simultaneously. (In fact, if solving just the LPA, rj K = 0; nevertheless, 
we have kept r\ K in eqs. (T2"2l - (|2"3"I) for a later use of these equations). When going to the second 
order, eq. f[2"3"|) acquires a new term and a new flow equation, that for p^, appears; and so 
on. According to Figure [fl an apparent convergence shows up. In all cases one observes 




FIG. 1: First four dimensionless fixed point couplings at zero momenta and zero external field: 
results obtained by truncating the flow equation, as a function of the order; the corresponding 
value for the complete equation is represented by the dotted-line. 

that: 1) there seems to be an oscillating convergence, 2) the value of /i w is found with about 
1% error truncating at order i + 3. 

Let us now turn to the study of the flow equation for the 2-point function coming from 

I II (2) 

the scheme proposed in [16(. As the effective potential (or the effective mass), r« p) can 
also be expanded in powers of the external field: 

00 cyn 

If P) = J2 TTTi Ti * n+2) (Pi °> °> • • • > °5 P)\r* ( 24 ) 

n=0 ^ U) - 

because 

r^(p,^0,0.-.-,0;p) = g|^ (25) 
10 




FIG. 2: Comparison of the self-energy when expanding only the flow equations for the self-energy 
p) and its derivative Z K (p) (strategy I): truncation is made at first (double dotted-dashed), 
second (dotted-dashed), third (dashed) and fourth (dotted) order; the complete solution is given 
by the straight line. In the figure, u = 5.9210~ 4 A. 

and we used that, at zero field, all odd vertex functions vanish. Equation (I24j) makes clear 
the point stated above: once approximation ([6]) is performed, truncating the expansion 
in powers of the external field is equivalent to neglecting high order vertices. Moreover, 
eqs. f l24l) and ( |25i) show that the procedure proposed in [l(| indeed includes all vertices, 
although approximately. 

We have now all the ingredients to discuss the main goal of this paper: the analysis of 
the expansion of the three flow equations for p K (p), Z K (p) and E K (p; p), eqs. (I17ti20p . around 
p = 0. In doing so, one can write: 

oo 1 
n=0 

and 

oo _ 

x^(p) = J2^x ( : ) P n . (27) 

n=0 

together with eq. (pill) . For example, when going to the first order, the six equations that 

11 



have to be solved are: 



and 



^4°>(p) = - (1 -;- /5) g (p) (28) 

(1 + H K ) 



K9«S^ ; (P) = (l+^)S^(p) + — ^y— / K (p,0) - — ^y— 



which correspond to the expansion of eq. (1201) . 

X 1 ) 



o (0) (0) (1 ~~ Vk/^Xk /Qn x 

= ( 30 ) 

and 



which correspond to the expansion of eq. ffTB]) . together with eqs. (1221) and (123]) . 

In fact, it is possible to perform two kinds of expansion. First, in order to isolate the 
effect of the field expansion just in the flow equations provided by the scheme presented in 
[if]], we shall expand only the flow equations for X K (p; p) and its derivative at zero momenta 
Z K (p), eqs. (1201) and ([TBI) , solving exactly the differential flow equation for p K (p), eq. (p3 
For example, at first order, one should solve simultaneously eqs. (I28H29I) . (I30H31I) . and (|Tj 
This is called "strategy I" . Second, to consider all the effects, we shall make the expansion 
in the three flow equations. For example, at first order, one should solve simultaneously 
eqs. fl2SH22P, fl3DffiHD, and (EMS}- We call this "strategy II". Notice that, as explained in 
171 ] . in order to get the correct scaling behavior it is mandatory to treat the equations for 
Z K (p) and S re (p;p) with the same approximations; it is then not possible to solve one of 
them completely while expanding the other one. 



Figure [2] presents the self-energy one gets 
I; it is also shown the function obtained in 



truncating up to fourth order, following strategy 



171 ] (from now on, the latter function, obtained 



by solving the 3 differential equations, eqs. (]T7]) . (jTSj) and f[2"U]) . shall be called the "complete 
solution"). Figure [3] presents the same results when following strategy II. These Figures 
show that, in both strategies of expansion, by truncating at first order one already gets a 
function with the correct shape in all momenta regimes. 
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FIG. 3: Comparison of the self-energy when expanding the three flow equations (strategy II): 
truncation is made at first (double dotted-dashed), second (dotted-dashed), third (dashed) and 
fourth (dotted) order; the complete solution is given by the straight line. In the figure, u = 
5.9210" 4 A. 

In order to make a quantitative evaluation of the approximate solution obtained doing 
the expansion, we have calculated different numbers describing the physical properties of 
the self-energy. First, as can be seen in both figures above, all solutions have, in the infrared 
(p <C u), the potential behavior characterizing the scaling regime: + p 2 ~ p 2 ~ v , where 
7] is the anomalous dimension. We have checked that, at each order and in both strategies, 
the resulting self-energy does have scaling, and we extracted the corresponding value of 
r). In fact, this can be done in two different ways: either using the ^-dependence of Z K 
(rj = — lim K _^ K d K log Z K ) or the p-dependence of S(p) stated above. We checked that those 
two values always coincide, within numerical uncertainties. Figure H] presents the relative 
error for rj, at each order, when compared with the value following from the complete 
solution. One observes: 1) in both strategies of expansion there is an apparent convergence, 
which is oscillatory; 2) the solution from strategy I reaches faster the correct result; 3) when 
following strategy I, already with a second order truncation the error is about 3% and it 
drops to less that 1% at the third order. Nevertheless, due to the mixed characteristic of 
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FIG. 4: Relative error (measured in percent) for the anomalous dimension, with respect to the 
value coming from the complete solution, as a function of the truncation order. Full line: expanding 
only the flow equations for the self-energy T, K (p; p) and its derivative Z K (p) (strategy I); dashed 
line: expanding all flow equations (strategy II). 

strategy I, when using this strategy at high order numerical problems arise: indeed, this 
task demands the numerical evaluation of high order derivatives of p K {p), to be used in the 
various flow equations obtained when expanding that of S K (p;p). If high precision in the 
result is required, strategy II is then numerically preferable. 

It is important to observe here that the procedure which can be compared to the usual 
truncation including a finite number of vertices is strategy II. Moreover, the inclusion of 
high order vertices without performing any other approximation is difficult; for example, 
the complete inclusion of the 6-point vertex has never been done. Accordingly, as can be 
seen in fig. HI when including only up to the 4-point vertex, as it usually done, the error in 
i] can be as large as 60%. 

A second number to assess the quality of the approximate solution is the critical exponent 
v. In order to calculate it, we extract the renormalized dimensionful mass from m 2 R = 
K 2 /i K (p = 0) and we relate it to the microscopic one by 

m 2 R {K = 0) ex {m 2 R {K = A) - m\ crit (K = A)) 2 ^, (32) 

where m^ cr u is the critical renormalized mass. With the complete solution one gets v = 
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FIG. 5: Relative error (measured in percent) for the critical exponent v, with respect to the value 
coming from the complete solution, as a function of the truncation order. Full line: expanding 
only the flow equations for the self-energy S re (p; p) and its derivative Z K (p) (strategy I); dashed 
line: expanding all flow equations (strategy II). 

0.647, to be compared with the best accepted value [31 ]: v = 0.6304 ± 0.0013. Figure [5] 
presents the relative error of the value of v extracted from the expansion. Once again, one 
observes that the convergence is much faster when following strategy I, i.e., when considering 
the effect of the expansion only on the self-energy equation. 

The large momenta regime (p u) of the self-energy can be calculated using perturbation 
theory, yielding the well known logarithmic shape: S(g) ~ Alog(p/B), where A and B 
are constants. For the complete solution presented in 17j one can prove analytically that 
A = m 2 /97t 4 , which is only 8% away from the exact result A = u 2 /96tt 2 (please observe that 
this coefficient is given by a 2-loop diagram for the self energy which is only approximatively 
included at this order). The proof of this analytical result remains valid when performing 
the field expansion, at any order and within both strategies. We have checked that our 
numerical solution always has the correct shape, with A = m 2 /97t 4 . This is due to the fact 
that already the first order in the expansion of E re (p; p) around p = contains the same 
2-loop diagrams contributing to the complete solution. 

In order to study the quality of the self-energy in the intermediate momenta regime, we 
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FIG. 6: Relative error (measured in percent) for A(cj) 2 ), with respect to the value coming from 
the complete solution, as a function of the truncation order. Full line: expanding only the flow 
equations for the self-energy T, K (p;p) and its derivative Z K (p) (strategy I); dashed line: expanding 
all flow equations (strategy II). 

have calculated a quantity which is very sensitive to this cross-over region: 



(the integrand is non zero only in the region 10 3 ^ p/u ^ 10, see for example [15(1 )• This 



quantity received recently much attention because it has been shown [32] that for a scalar 
model with O(N) symmetry, in d = 3 and N = 2, it determines the shift of the critical 
temperature of the weakly repulsive Bose gas. It has then been widely evaluated by many 
methods, for different values of N, in particular, for N — 1. With the numerical solution 

n 

found in [17] . one gets a number almost within the error bars of the best accepted results 
available in the literature, using lattice and 7 loops resumed perturbative calculations. Please 
observe that these errors are as large as 10%, which is an indication that this quantity is 
particularly difficult to calculate. In Figure [6] we plot the relative error in A(0 2 ), at each 



order of the expansion, when compared with the complete solution result found in [17J. One 
can appreciate that 1) for both expansion strategies there is an apparent convergence, which 
is also oscillatory; 2) in both strategies, already with a second order truncation the error is 
about 1%. 
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IV. SUMMARY AND CONCLUSIONS 



In this article, the inclusion of a finite number of vertices in NPRG flow equations is 
analyzed. An unsolved difficulty of this usual strategy (originally proposed by Weinberg) 
is the estimation of the error introduced at a given step. Moreover, without performing 
further approximations, it is very hard to reach high orders of the procedure. The study of 
its convergence is thus a difficult task. In the present work we analyse this problem using a 
different approximation scheme [if]]: instead of considering a finite number of vertices, this 
procedure includes all of them, although approximately. Within this context, it is possible 
to estimate the error of the Weinberg approximation, order by order. To do so one can 
perform, on top of the approximation presented in [3] , the usual truncation in the number 

the 2-point function of 



171 ] that, at least in this 



of vertices. The analysis has been done in the particular case o 
the scalar field theory in d = 3 at criticality. It has been shown 

case, the procedure proposed in Q yields very precise results. Auother iuterestiug outcome 
of the present work follows from the fact that, within the approximation [la], truncation 
in the number of vertices is equivalent to an expansion in powers of a constant external 
field. The latter is usually employed in the DE context in order to deal with complicated 
situations. The analysis of the present paper generalizes this expansion procedure when non 
zero external momentum are involved. 

The calculation of the 2-point function demands the study of both the p = and the 
p 7^ sectors. While the first one is given by the well studied DE flow equations, the latter 
follows from the approximation scheme introduced in 16| to calculate the flow of p). 
We used two different strategies to perform the field expansion, both of them around zero 
external field: either expanding only the flow equation for the self-energy (and its derivative) 
(strategy I), or both the effective potential and the self-energy (and its derivative) flow 
equations (strategy II). We have studied the convergence of various quantities measuring 
physical properties of the self-energy in all momenta regimes: the critical exponents rj and 
v of the infrared regime, the coefficient of the ultraviolet logarithm, and A(0 2 ) which is 
dominated by the crossover momenta regime. 

As stated in section 1111} the strategy that can be compared to the usual truncation which 
includes a finite number of vertices is strategy II. For example, including completely the 
4-point vertex as it is usually done (i.e., in the language of field expansion, going only up to 
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the first order of the expansion), when describing the deep infrared regime one could make 
errors as big as 60% in the critical exponent 77 (see figure H]). If one wants results with less 
that 5% error for this quantity, the inclusion of up to 8-point vertices (i.e., going up to third 
order) is necessary. 

However, when going to higher orders in the field expansion, the series for all considered 
quantities seem to converge rapidly, within both strategies. The convergence is faster when 
using strategy I, i.e., when making the expansion only for the approximate flow equation 
resulting from the method presented in 16j. For example, using strategy I, a third order 
truncation introduces a relative error smaller that 1% for all studied quantities; while using 
strategy II, in order to reach the same error one needs 6th order for 7], 4th order for v and 
2nd order for A(<f) 2 ). Nevertheless, due to numerical difficulties, if trying to go to high order 
expansions, it is preferable to use strategy II, i.e., expanding also the effective potential flow 
equation. 

It is difficult to assess the generality of these results on the use of field expansion on top of 
the strategy proposed in 16J]. Of course, there are situations where expanding in an external 
field is not a priori convenient. One can mention as a first example, situations where there 
is a physical external field (as in a broken phase or when an external source for the field is 
considered). A second example is two-dimensional systems where even in the DE, the field 
expansion does not seem to converge. Nevertheless, the short study presented in the present 
paper allows to consider field expansion on top of the approximation proposed in [16J as a 
possible strategy to deal with many involved models, as for example QCD. 
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